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Abstract 

We apply the finite-element lattice equations of motion for quantum elec- 



(N 
> 

in 

C^ . trodynamics given in the first paper in this series to examine anomalies in 

. the current operators. By taking explicit lattice divergences of the vector and 

Qv . axial-vector currents we compute the vector and axial- vector anomalies in two 



and four dimensions. We examine anomalous commutators of the currents to 



Oh. compute divergent and finite Schwinger terms. And, using free lattice propa- 



gators, we compute the vacuum polarization in two dimensions and hence the 



^ , anomaly in the Schwinger model. A discussion of our choice of gauge-invariant 

current is provided. 
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I. INTRODUCTION 

This paper summarizes the status of the finite-element approach to gauge theories, in 
which the Heisenberg operator equations of motion are converted to operator difference 
equations consistent with unitarity. For a review of the entire program, from quantum 
mechanics to quantum field theory, see Ref. |T|. The gauging of an Abelian theory on a 
finite-element lattice was first carried out in Ref. [^, while the generalization to a non- 
Abelian theory was given in Ref. [Q . The beginnings of constructing a canonical quantum 
electrodynamics along these lines were given in Ref. [Q. In particular, we showed there that, 
at least for background fields, the Dirac equation is unitary, in the sense that canonical 
anticommutation relations hold at each lattice time. 

In this paper, we will first show, in Sec. ||, how to construct free lattice propagators 
for electrons and photons. We remind the reader how interactions are introduced, in Sec. 



Ill, and then, using the gauge- invariant current, compute directly vector and axial- vector 



anomalies. Another anomalous aspect of the current operators is examined in Sec. |^, 
where equal-time current commutators are computed, from which the divergent and finite 
Schwinger terms are computed. In Sec. a lattice loop calculation is described. We compute 
the vacuum polarization in two dimensions using the electron propagator found in Sec. 0, 
and obtain results consistent with the known continuum Schwinger model. The continuing 



direction of our program is sketched in the concluding Sec. 0. In the Appendix we discuss 
a possible alternative definition of the current possible in Euclidean space, and show why 
such a current is unacceptable. 



II. FREE DIRAC EQUATION. LATTICE PROPAGATORS 

We begin by reminding the reader of the form of the free finite-element lattice Dirac 
equation^: 

Here /x is the electron mass, h is the temporal lattice spacing, A is the spatial lattice 
spacing, m represents a spatial lattice coordinate, n a temporal coordinate, and overbars 
signify forward averaging: 

By transforming this into momentum space we find that the momentum expansion of the 
canonical Dirac field is 



y u! 



(2.3) 



where the spinors are normalized according to 



ES±4/ = T^-±A., (2.4) 

s 

where u = u_ and f = «+, and where 

2t 
p = — , u; = / = v/^^TT^, (tp), = tanpivr/M. (2.5) 

In (|2.3|) A is the eigenvalue of the Dirac transfer matrix 

A = i±^ ^ e^-('^)^ (2.6) 

1 - ihuj/2 ^ ' 

Notice that we may solve ( p.6|) for lo: 



^Note that we use a different sign of ii than in Ref. |^, and correspondingly, a different definition 
of spinors in (|2.4|) and a change of sign in the unitary time evolution operator [/, (|3.2|). 



2 hn 

ijj = —tan — . 
h 2 



(2.7) 



We take M, the number of lattice points in a given spatial direction, to be odd, so that ip is 
periodic on the spatial lattice. It follows from ( p.3|) that the canonical lattice anticommuta- 
tion relations for the Dirac fields 



{V'm,n'^m',n} = ^'^m,m' 



(2.8) 



are satisfied if (L = MA, the length of the spatial lattice) 



L3 



i"p,S5 "p'.s'l ~ TS^P^P' 



Ss,s', 



(2.9a) 
(2.9b) 



and all other anticommutators of these operators vanish. 

At this point we can construct a free lattice Dirac Green's function according to 



(2.10) 



where the step function is defined by0 



ri{x) 



1, X > 0, 
1/2, x = 0, 
0, x<0. 



(2.11) 



Using the fact that b and d are annihilation operators, we easily see from the Fourier expan- 
sion ( p.3|) that the free Green's function has the following momentum expansion: 



^Inclusion of the x = value of the step function eliminates the following nonlocal term from 

O): 

^n,n' :^ E e^P-('"-'"')2-/^(/i - 7 • P)/ tan hi}/2. 



P 



+ A_A"-"'e-'P('"-'"')2-/^r7(n' - n)]. (2.12) 

To recast ( |2.12| ) in four-dimensional momentum space, we can use the identity for integer p 



i—- = J ' ^- ' (2.13) 

"'T 0, p < 0, 

where 7 is a contour encirchng and 1 in a positive sense. It is then a straightforward 
exercise to show that, apart from a contact termQ, 

'j'^ sin hfl + [fi — J ■ p)hcos^ hCl/ 2 
COS h{Q — ie) — cos hO, 

This expression has the expected continuum hmit: as /i ^ and nh —>■ t, ( |2.14|) becomes 



— r rff]e-*^(*-*')^ ^e^P(— ')2-/M_I:l_Ji_, (2.15) 

where p = (fi, p). 

In a similar way we can work out the free lattice Green's function for the photon. We 
first note, from the lattice versions of the Maxwell equations in the temporal gauge, Aq = 0, 
in part, 

E = A, E + VxB = 0, (2.16) 



namely. 



1 



^m,n — ,\-^m,n+l -^m,n)j [^■^') 



^This contact term is 



^ <5„,n'E^^^'^'""'"'^'^^''(^-^-P)' 



4L3 

p 



and serves to convert, in the numerator, cos hil./2 to cos /ifl/2 



and 

T[{Ei)m,n+l - (-^j)m,n] + -^^ijk[{Bk)mj+l,m^,n " {Bk)mj,m^,n] = 0, (2-18) 

that the momentum-space transfer matrix for the photon is 

T = e'^^ (2.19) 

where 

2t 2 hn , , 

just as in ( |2.7| ). [Incidentally, notice that this form of the free photon transfer matrix 
establishes unitarity of the free lattice Maxwell equations.] 
Now we compute the photon propagator defined by 

= {Al^^^AL^,Jr^{n - n') + (A^, „,Ai,,„)r^(n' - n), (2.21) 

which, on use of the canonical commutation relations for the photon creation and annihi- 
lation operators, defined through the momentum expansion of the photon fields, (2.17) of 
Ref. i, 

K, 4,] = 0, K, 4!] = 5k,k'r(k), (2.22) 

where the transverse projection operator is 

r (k) = 5., - ^^j^, (2.23) 

leads to 

L^ ^ 2ujk 

k 

Use of the identity ( p.l3|) leads to the four- dimensional form, apart from a contact term. 



6 



/)2 r-^/h , 1 



ik-(m— m')27r/M eij 



r(k) 



X- 



cos^ /ifi/2 



cos hiyt — ie) — cos hVt 
completely analogous to the electron Green's function ( |2.14| ). 



(2.25) 



III. INTERACTIONS. AXIAL- VECTOR ANOMALY 



Interactions of an electron with a background electromagnetic field is given in terms of 
a transfer matrix T: 



^/'n+l = Tni^n, 



(3.1) 



which is to be understood as a matrix equation in m. Explicitly, in the gauge A*^ = 0, |^ 
T = 2f/-i - 1, f/ = 1 + "^ - A/^ . X>, (3.2) 



A 



where 



^m m' = -(-1)"'^+"'^ [e^,,m'. COS Cm,,m'. " ^ slu Cra„m'] SBC (^^^ 6, 



Here 

1, -m > m', 
Cmm/ = N 0, -m = m', 
—1, vn < m', 

and (the following are local and unaveraged in m±, n) 

eA 



J' ji 






ruj, 



and 



A/ 



Cm,,mr = ^ sgn (m^' - mj)sgn {m'- - m'j)C„ 



S'=i 



(3.3) 



(3.4) 



(3.5a) 



(3.5b) 



with 

sgn (m - m') = em,m' - Sm,m'- (3.6) 

Because T> is anti-Hermitian, it follows that T is unitary, that is, that <^m,n = V'm,n is the 
canonical field variable satisfying the canonical anticommutation relations (|2.8| ). 



It is instructive, and very simple, to consider the Schwinger model, that is the case with 
dimension d = 2 and mass fi = 0. Because the light-cone aligns with the lattice in that case, 
we set /i = A. Then we see that the transfer matrix for positive or negative chirality, that 
is, eigenvalue of ^75 = •y^'j^ equal to ±1, is 

r. = j^. (3.7) 

From (|3.3| ) we see that the numerator of T is 

(1 ± P)„,„, = \bm,^.e^'^ T (-l)"+'"'e,n,,n'e-*^'"^'"'^™-'] sec C, (3.8) 

while it is a simple calculation to verify that the inverse of this operator is 

(1 + ^);;!,n' = \ (^m,-' + '^m,m'-ie-'''^-') , (3.9a) 

(1 - ^);^Im' = \ ('^m,™' + 5™,™'+ie'^^-) . (3.9b) 

It is therefore immediate to find 

(-'+)m,m' = Orra,m'+lC " ; (3.10a) 

(T_)™,^, = (5^+i,^,e-2<'"', (3.10b) 

which simply says that the + (— ) chirality fermions move on the light-cone to the right (left), 
acquiring a phase proportional to the vector potential. Solution ( |3.10| ) directly implies the 
chiral anomaly in the Schwinger model, as shown in 0. 

To compute the vector and axial-vector anomalies in (3 + 1) dimensions, we will find it 
useful to expand T in powers of the temporal lattice spacing h (but not in the spatial lattice 
spacing A). This leads to a very simple form for the transfer matrix: 

8 



T = l + h 



-^f,JO ^ _^o^ . ^ 



h' 






0(/i3). (3.11) 



On the lattice the gauge-invariant current is written in terms of the all-averaged field \&, 

1 



mi" Q VV^m,"+l ' r'm,nj; 



rather than the canonical field 6. That is, the current is 



JiL,n = e^U7V^^,n = 4[0Ul + Ti)]„.7V[(l + TJ0J, 



(3.12) 



(3.13) 



(For a discussion of why this choice of current is used, see the Appendix.) In the absence of 
interactions it is easy to show that 



("9^JM") = 0, 



(3.14) 



where the quotation marks signify a finite-element lattice divergence. At the initial time 
n = we introduce the momentum expansion ( ^.3[ ), and compute the expectation value in 
the corresponding Fock-space vacuum. If we use the expression for the current in ( ^.131 ) and 
keep only terms of 0{h), we find 

(^^,i) = (Jko) - E x(<^L'.o7°7 ■ kl^m + ^(l))^^^7V0„.,o) + h-c. (3.15) 



Thus, to leading order in h, the vacuum expectation value of the divergence of the current 
is (averaging in m is implicit) 



':%j''n 



h 



2e 



AL3 



5Z 5Z ~ sin[pj(mj. - mj)27i/M] 



rrij ,m ■ ' 



(3.16) 

m' p 

where a = — ImP. Not only is this not zero, it is not gauge invariant, for under a gauge 
transformation, 



6V' , 

rrij ,m - 



ie{6Qrn, -6Qrn',)1^i 



<■]'"•'] 



(3.17) 



This seems to be a conundrum, because the current was explicitly constructed to be gauge 
invariant. The resolution lies in the recognition that since the canonical field is covariant 
under time-independent gauge transformations. 



^0m,n = ieSQm4>in,m (3.18) 

the states created by that operator must transform hkewise. That is, we require 

H<I^UK'} = *e(5fi™; - (5fi„J(0j„^,0^,). (3.19) 

This suggests we supply the following factor in computing the vacuum expectation values: 

because 

OC,mi,m' ^ C(,Oi'm,, Oi 6^' j^m' ,m,- ■ 1,3. ^ii) 

■'J ■' J J -^ 

Then, the derivative operator that occurs in (|3.16| ) becomes 



-Ime'^'"^''"^''"^''"^P^' ,=0, (3.22) 

since the term with ttlj = m', does not contribute. 

In the next order in h it appears that a nontrivial contribution could emerge. It is the 
following gauge-invariant structure (here only the the leading term in AA is written down): 

("V'^")(x-^M^C«'- (3.23) 

i 

This is a nonlocal term without a continuum analogue, because its gauge invariance depends 
on the periodic boundary conditions on the spatial lattice. However, to this order, since 
^m,i = ^m,o + 0{h), a simple inspection shows that the coefficient of this term is zero. 



The axial current is defined just as in (|3.13| ) with the additional factor of 275. We compute 



the divergence analogously to (|3.16| ), except here the first apparently nonzero term is 0{h^): 



oc--^ y y^e,,,sin[p,(m'-m'')27r/M]«(') , a^'^ „. (3.24) 



Because of the diagonal property of V in the perpendicular coordinates [see (PTB|)], ( |3.24| ) is 



zero. So, to 0{h) both the vector and axial- vector currents are nonanomalous. This may be 
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a bit surprising since when h = A a.n axial- vector anomaly emerges in (1 + 1) dimensions. We 
are presently trying to understand this state of affairs by carrying out an exact calculation 
with /i = A in (3 + 1) dimensions. 

IV. ANOMALOUS CURRENT COMMUTATORS 

It is extremely interesting to compute commutators of the current defined by (|3.13|) , and 
compare with the anomalous commutators in the continuum PJ^: 

([J°(0, x), J(0)]) = tSV6{^) + tdVV^6{^), (4.1) 

where S is the quadratically divergent Schwinger term, and d = l/127r^ for the Bjorken- 
Johnson-Low regularization. Straightfoward evaluation of the current-current commutators 
on the lattice requires use of the transfer matrix T for the free Dirac field: 

'/lA 2j \ h A2 

1 + V + ^t^ J (l - V - A^' + A ^^ ^ • * " ^^'^) ' ^^-2) 

where t = tp is given by ( [^.5| ). The result is the following expression for a hypercubic lattice, 

h = A: 

njo p n 4e^ V tan(7rA:;./M) 2.^k.(m-m')/M 

\[ m,n, m',nj; ^^ ^y Z. (1 + ( A^k-k'/2)2) (1 + (Aa;k'/2)2) (AcUk'/2) 

(4.3) 

This matrix (in m and m') is shown for lattice size M = 9 and mass /i = in Fig. |l|. The 
abscissas of the plot correspond to one-dimensional representations of the three components 
of m and m', respectively. (That is, the base-M number (mi, m2, m^) is converted to a base- 
10 value of the abscissa.) It is apparent that the leading-order behavior of this commutator 
is a first derivative of the delta function, as expected. 

To make this comparison quantitative, we fit this result to the functional form expected 
in the continuum. The lattice analog of the Dirac delta function is 

11 



so we take as a trial function 

When /i = 0, the coefficients a do not depend on the lattice spacing A; they will presumably 
remain finite in the continuum limit, A — > 0. The functions Pr(^) are inserted to force each 
term in the series to be of the same order in A as the commutator ( |4.3| ), namely 1/A^: 

P,(A) = A^'^-l (4.6) 

Therefore, in the continuum limit, the first term in the series (|4.5| ) will be quadratically 
divergent, the second will be finite, and the rest will vanish, as expected. 

To perform the fit, the problem is first converted from a three-dimensional to a two- 
dimensional representation. (This is necessary because the generation of three-dimensional 
plots is impractical beyond M = 9 and three-dimensional fits are prone to difficulties.) 
Specifically, a plot such as Fig. |l| is projected onto a plane orthogonal to the m, m' plane 
and to the diagonal of the matrix. The resulting curve is then fit to a similar projection of 
( [4.5|) . In Fig. H we show the results for the coefficients of a fit including only the first three 
odd derivatives of the lattice delta function (that is, for i? = 2 in (|4.5| )) for various lattice 
sizes. The coefficient of the first derivative is roughly consistent with the Schwinger result 
0. Similarly, the result for the coefficient of the third derivative term in the commutator is 
in approximate agreement with the agreement with the BJL result 0; but seems to be closer 
to that of Chanowitz |]^. It is not surprising that the results are regularization dependent, 
given the divergent nature of the commutator. A more detailed discussion of this calculation 
will appear elsewhere 0. 

V. LOOP CALCULATION. SCHWINGER MODEL 

It is illuminating to calculate physical quantities such as the magnetic moment or the 
axial-vector anomaly by doing perturbation theory on the lattice. That is, we calculate 

12 



quantities such as vacuum polarization by using weak-coupling perturbation theory together 
with the lattice propagators given in ( |2.14| ) and ( |2.25| ). We illustrate this idea in the simple 
context of the Schwinger model, two-dimensional QED. (An analogous calculation in four 
dimensions of both the axial-vector anomaly and the anomalous magnetic moment is in 
progress.) 

We begin by discussing the Schwinger model {d = 2 electrodynamics) in the continuum. 
To calculate the anomaly it is essential to work with massive QED2 and then take the 
massless limit to avoid infrared singularities. The vacuum polarization is then 



We require that 11'^'" be transverse: 

= nrPVr, - - 

A;2 



W^ = U{e){g,^-^). (5.2) 



Contracting on the indices gives 

n(A;2) = -^e^ J (27,)2 ((^ _ A:)2 + ^2)(p2 + ^2)' (5-3) 

which is easily expressed using Feynman parameters as 

n(^') = — / dx- '^ ^ . (5.4) 

For the /x 7^ theory, we require that 11(0) = 0. This implies the presence of a subtraction 
constant e^/vr, which for the fi = theory is the square of the boson mass, or the anomaly. 
Rather than introduce a Feynman parameter we could obtain this result by setting k = 
in ( p.3| ) and performing the momentum integration explicitly. We have, then, for the 
unsubtracted vacuum polarization 

e^ f 4/i2 eV^ /"~ 1 
U{k = 0) = -ijT^^ / dpodpi-^ 2 2 I M2 = / dpi-- -y, (5.5) 

(27r)2y (p2 _ p2 _ ^2 + ^g)2 ^ J^ (p2^^2)5 

which again gives the usual result, 11(0) = — c^/tt. 

We now describe the analogous calculation in the finite-element lattice framework. The 
lowest-order vacuum polarization is given by the time-ordered products of currents, 

13 



where the gauge invariant electromagnetic current is given by ( |3.13| ). Using the lattice 
Green's function ( ^.141 ) and the transfer matrix ( [4 .21 ), we find 

^^-ip,im-m')2./M^ip,{m-n.')2^/M^^^^I.^^ (^7° sin(/ini) + (/i - -iPi)hcOs''{h(li/2)\ 



(5.7) 



xPi7"7^2 (^7 sin(/ifi2) + (/i - 7P2)/icos'(/ifi2/2)j P2] 

1 

X — , 

{cosh{Vti — it) — coshVti){cosh{Vt2 — it) — coshVt2) 

where P = (1 + T)/2 and V = 7''P^7°. We now perform a lattice Fourier transform: 

M 00 

LU 2 



n'^^(g, .;) = A E E ^ ^i^ ^e-^^-("-"')e-(™-')^-/^C^„^^, „,. (5.8) 

m=l n=— 00 



Note that 

M 



E e^«™2"/^^ = M(5g,o and E e"'''"" = 27r5(/icj), -tt < /icj < tt. (5.9) 



m=l n=— 00 

Upon contracting with g^j^^ we find for the vacuum polarization 

n g, u; = -ie^ - — 27rAM ^ sm — 
^ /"^"/^ rff]2l6z/Vcos2(/ifii/2)cos2(/ifi2/2) 



_ (cos/i(fii — it) — cos/irii)(cos/i(i72 — it) — cos/ir22) 
1 
^ (1 + rV2 + r2t2)(i + r2z/2 + rH^) ' ^^'^^^ 

where u = |/iA and r = h/A and where, in the trace, we recall that 7'^7"7a = in two 
dimensions. Here we understand that Qi = uj + CI2 and ti = tp^ = tg+pj. The integral over 
CI2 can now be done to yield, when tg = uj = 0, 

u(nn\ = -^\^A r' ly' cos' hn2/2 cos/i^2 / 4sin^/i^2/2 \ 

^HU,Uj vr ^ M(l + r2i.2 + r2t2)2sin3/,fi2 V cos/ifis J' ^ ^ 

For zero external momentum, fii = f22, where 

/ifi2 = 2tan-^(r(t^ + z/2)5). (5.12) 
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We perform the momentum sum in (|5.11|) for various values of the lattice size M, mass 
parameter z/, and ratio of temporal and spatial lattice spacing r. The parameter v is con- 
strained by M~^ ^ z/ ^ 1. Figure ^ is a graph of the sum for M = 2533 versus z/. Plotted 
are curves with r = 0.1, 0.75, 0.99. The first two cases stay close to 1 for a significant range 
of u, and only for r = 0.99 does the value drop significantly below 1. Note that for very small 
u the sum diverges, as one expects for u ~ M^^. This feature is shown more explicitly in 
Fig. |, that is a plot for r = 0.1 for lattice sizes 1105 < M < 3145, with increasing endpoint 
for decreasing lattice size. For smaller lattices, the curve turns up for larger i>, as the limit 
of the domain of validity of the lattice calculation is reached. 

More detailed analysis shows that, as a function of r, —11(0, 0)/(e^/7r) is on the average 
unity, but is very noisy, with wild fluctuations around 1. For small r, fluctuations are greatly 
diminished, but for r ~ 1 and for various values of the input parameters, hQ2 can get very 
close to vr and make a few terms in the sum very large. (One might think that the desired 
value of r is r ~ 1, the case of a hypercubic lattice; as noted in Sec. |T| the anomaly in the 
Schwinger model has previously been calculated with the finite element lattice technique 
with r = 1.) This question of stability is also under intense scrutiny. 

VI. CONCLUSIONS 

This paper is a report of recent progress in our finite-element program as applied to elec- 
trodynamics. More complete treatments of the analyses of the last two sections will appear 
elsewhere |P,0. Our immediate goal in QED remains the extraction of a nonperturbative 
value for the anomalous moment of the electron. 

Of course, our real interest lies in non-Abelian theories. Our first task there is to complete 
the formulation in four dimensions PJTO|]. This should be accomplished in the next few 
months. Then we can apply this approach to theories such as QCD. 
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APPENDIX: CURRENT CONSTRUCTED FROM EUCLIDEAN LAGRANGIAN 



In the text we simply assumed a form of the current (|3.13|) which was manifestly gauge 
invariant. We are at liberty to do so, because the Minkowski finite-element equations of 
motion are not derivable from a Lagrangian. The current cannot be derived from the Dirac 
equation, but is an independent source for the Maxwell equations, see Ref. [^. However, if 
one were to work in Euclidean space-time (periodic in all four directions), it is possible to 
construct an action from which the equations of motion are derivable, and which therefore 
supplies a lattice current. The fermion part of that action is (a factor of i is absorbed in 
going to Euclidean space) 



iy/ = A^ J]^;,(^|r.x» + /.) 



^^„ r = (7V.,7°). (Ai) 



0,, „,o 

m'- ^ = ''' 
m,m' 



If we vary ( |XT| ) with respect to ip"^ we obtain the Euclidean Dirac equation 

|-r.X> + /i')^ = 0. (A2) 

where, as in (|Al|), a four-dimensional scalar product is implied. 

Given an action, we can construct a conserved vector current by making a local gauge 
transformation 

5^jn = ie5n^ijjjn- (A3) 

Because the Dirac equation, and hence the action, is invariant under the global version of 



(|A3|), 6^ni = S^ = constant, we must have by the action principle 
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5Wf = Q = -Y, Jln\mrn.,m^ - 5fi™.-l,m J, (A4) 



m 

from which we read off the conserved current 



[mi - m") 



JL = —e y ih—, — V^ih—i, — sernfmj — m') sm( 

m'- ,m'' 

x(-l)™^+<secC(*)exp(-ze™;.W'Cm'<). (A5) 



z' i -' z' z 



(The same result, of course, can be obtained by varying V with respect to A ^ _-^ ^ .) The 
expression for this Euchdean current has been simphfied by deleting constant terms. It is 
easy to verify explicitly that this current is both conserved and gauge invariant. Similarly, 
by making a chiral transformation, 

Si^m = 7°75^^m^m, (A6) 



we can construct the axial- vector current J^^, which has the form of ( |A5| ) with the replace- 
ment 

By construction, these currents possess no anomalies. However, they appear to be completely 
unacceptable, because they are horribly nonlocal. In particular, they possess no Minkowski 
analogues, in the sense that it is not possible to analytically continue back to real unbounded 
times. Crucial to our formulation is the propagation of the operators from past times to 
the present time, so that we can solve for the field operators by time-stepping through the 
lattice. The Euclidean current (^) involves fermion field operators at all Euclidean times. 



which would make it impossible to solve for the operators at time n in terms of operators 
at earlier times. Therefore, for the considerations of the text we use the gauge-invariant 
current (|3.13|) and its axial analogue, currents which do possess anomalies. 
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FIGURES 
FIG. 1. Three-dimensional plot of ([4 .31) for M = 9 and /x = 0. 



FIG. 2. Coefficients of three spectral fit for different lattice sizes (/x = 0). For the first derivative 
term, the coefficient shown is S'A^. 



FIG. 3. Plot of the sum (5.11) for M = 2533 as a function of z^ = /iA/2. Shown are curves 



with r = h/A = 0.1, 0.75, 0.99. 



FIG. 4. Plot of ( ISHP for 1105 <M< 3145 and r = 0.1. 
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